Transport Properties of Anisotropic Polar Fluids: 
1. Quadrupolar Interaction 

G.A. Fernandez, J. VrabecE], and H. Hasse 
Institute of Thermodynamics and Thermal Process Engineering, 
University of Stuttgart, D-70550 Stuttgart, Germany 



Number of pages: 45 
Number of tables: 1 
Number of figures: 14 



* To whom correspondence should be addressed, tel.: +49-711/685-66107, fax: +49- 
711/685-66140, email: vrabec@itt.uni-stuttgart.de 



ABSTRACT 



Equilibrium molecular dynamics simulation and the Green-Kubo formalism 
were used to calculate self-diffusion coefficient, shear viscosity, and thermal 
conductivity for 30 different quadrupolar two-center Lennard- Jones fluids along 
the bubble line and in the homogeneous liquid. It was systematically investi- 
gated how anisotropy, i.e. elongation, and quadrupole momentum influence the 
transport properties. The reduced elongation L* was varied from to 0.8 and 
the reduced squared quadrupole momentum Q*"^ from to 4, i.e. in the entire 
range in which parameters for real fluids are expected. The statistical uncer- 
tainty of the reported data varies with transport property, for self-diffusion 
coefficient data the error bars are typically lower than 3%, for shear viscosity 
and thermal conductivity they are about 8 and 12%, respectively. 

KEYWORDS: Green-Kubo; molecular dynamics; quadrupole; self-diffusion; 
shear viscosity; thermal conductivity. 
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1 INTRODUCTION 



From the pioneering work of Alder et al. |;1.2J who used molecular dynamics 
to investigate transport properties of hard-spheres, this method has proved to 
be a successful tool for the development and test of theories [3J ; examples are 
the discovery of the hydrodynamic long-time tail [4j or the development of 
advanced kinetic theories [SHU] . Due to the increasing computer power, molec- 
ular dynamics is today an interesting option for studying [7] and predicting 
[8] transport properties. In particular, for predictions of transport properties 
of liquids, where no satisfactory analytical theory exists, molecular dynamics 
is the most suitable available method. 

Molecular simulation is also attractive because it allows rigorous testing of 
theories and models, and systematically studying the influence of any molecu- 
lar parameter on any transport property. Although real fluids usually consist 
of polar non-spherical molecules, most extensive studies on transport proper- 
ties have been done on the basis of very simple molecular models, e.g. hard 
sphere potential PIU] , square well potential |10|lllll2lll3fl4|15j . steeply repul- 
sive potential (TB], or spherical Lennard- Jones potential P^fTTj . 

In the present work, a systematic study on the influence of elongation and 
quadrupole momentum on self- diffusion coefficient, shear viscosity, and ther- 
mal conductivity of two-center Lennard- Jones plus point quadrupole (2CLJQ) 
fluids is carried out in the liquid region covering a broad range of temperature 
and density. This work is based on the knowledge of accurate vapor-liquid 
equilibria of 2CLJQ fluids from previous publications of Stoll et al. [T8IIT9] . 
Molecular simulations are carried out along the bubble line and in the homo- 
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geneous liquid for different models over a grid of reduced temperatures and 
densities. Thus, a consistent comparison of the transport properties of different 
models is possible and subsequently the effect of elongation and quadrupole 
momentum can be identified. 

It has been recently shown that the 2CLJQ model is not only an interesting 
model fluid but also suited for accurately describing the properties of real 
fluids. Vapor-liquid equilibria of 35 pure substances were successfully mod- 
eled with that approach [20ll2T] . Also for mixtures good results were achieved 
P!^22f23j . Furthermore, properties such as Joule-Thomson inversion j23], self- 
diffusion and binary Maxwell- Stefan diffusion coefficients [25j, shear viscosity, 
and thermal conductivity [26] were reliably predicted by 2CLJQ models which 
were parameterized using vapor-liquid equilibrium data only. 



2 MOLECULAR MODEL 



The intermolecular interactions are represented by the two-center Lennard- 
Jones plus point quadrupole (2CLJQ) potential. The 2CLJQ potential is pair- 
wise additive and consists out of two identical Lennard- Jones sites a distance 
L apart (2CLJ) plus a point quadrupole of momentum Q placed in the geo- 
metric center of the molecule oriented along the molecular axis connecting the 
two Lennard- Jones (LJ) sites. The interaction energy of two molecules i and 
j is given by 

\rab/ \rab/ J 

Here, Tab is one of the four LJ site-site distances; a counts the two LJ sites of 



a=l b=^ 
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molecule i, b counts those of molecule j. The LJ parameters a and e represent 
size and energy, respectively. The contribution of a point quadrupole is given 

by m 



1 3 Q 
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4:7160 4: IVijl 



1 - 5 (c2 + c^) - ISc^c^ + 2 {siSjC - AciCjf , (2) 



with Cfc = cos9k, Sk = sinOk, and c = coscpij. Herein, is the center-center 

distance vector of two molecules i and j. 6i is the angle between the axis of 

the molecule i and the center-center connection line and (pij is the azimuthal 

angle between the axis of molecules i and j. Finally, eo is the electric constant 
in 8.854187817-10-12 Q2/p 

A specific 2CLJQ model, e.g. for a real fluid like nitrogen, is fully determined 
by five parameters: a, e, L, Q ^UfIT\ and the molecular mass m. But in 
molecular simulation all relevant physical properties can be treated in a re- 
duced form. Here, they are related to a, e, and m, so that the reduced results 
are valid for all combinations of these three parameters. In this form, only two 
molecular parameters remain, i.e. reduced elongation L* = L/a and reduced 
squared quadrupole momentum Q*'^ = Q'^ / (Aneoea^) . Henceforth, "squared" 
will be omitted in the text for brevity. 



3 TRANSPORT COEFFICIENTS 



Transport properties were calculated by equilibrium molecular dynamics sim- 
ulation and the Green-Kubo formalism pSf^ . In this formalism, transport 
coefficients are obtained by integrating time autocorrelation functions of the 
corresponding microscopic fluxes. 
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The self-diffusion coefficient D is a measure for the mobihty of individual 
molecules within a fluid. It is calculated by integration of the single molecule 
velocity autocorrelation function [30f3T] 

^^3nI (vfc(t)-Vfe(O)), (3) 

where Vfc(t) expresses the velocity vector of molecule k at some time t, and 
<...> denotes the ensemble average. Eq. yields the self-diffusion coefficient 
averaging over molecules. 

The shear viscosity 1], as defined in Newton's "law" of viscosity, describes the 
resistance of a fluid to shear forces. It refers to the resistance of an infinitesimal 
volume element to shear at constant volume [32] • The shear viscosity can also 
be related to momentum transport under the influence of velocity gradients. 
From a microscopic point of view, the shear viscosity can be calculated by 
integration of the time-autocorrelation function of the off-diagonal elements 
of the stress tensor, i.e. J^'^ [30|3T] 

where V is the molar volume, kB is the Boltzmann constant, T the tempera- 
ture, and < ... > denotes the ensemble average. The statistics of the ensemble 
average in Eq. (jlj) can be improved using all three independent off-diagonal 
elements of the stress tensor, i.e. J^^, J^^, and J^^. For a pure fluid, the 
component J^^ of the microscopic stress tensor Jp is given by pS] 

N 1 TV AT 3 3 Q„, 

J- = E - ^ E E E E -^P- (5) 

i=l ^ i=\ j+i k=l 1=1 ^' kl 
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Here, i and j are the molecular indices. Lower indices / and k count all sites, 
including the quadrupolar site, and the upper indices x and y denote the vector 
component, e.g. for velocity vf or center-center distance rfj. 

The thermal conductivity A, as defined in Fourier's "law" of heat conduction, 
characterizes the capability of a substance for molecular transport of energy 
driven by temperature gradients. It can be calculated by integration of the 
time-autocorrelation function of the elements of the microscopic heat fiow 
and is given by [30ll3T] 



V/f dt(j,^(t)-J3)). (6) 



The expression for the heat fiow 3g in pure fiuids has been derived by Evans 
[33] and is given by 

^ i=l j^i ^ i=l k=l 1=1 ^'^f'l 

(7) 

where Wj is the angular velocity vector of molecule i, Ij its matrix of angular 
momentum of inertia, and Uij the intermolecular potential energy. The torque 
Fij refers to a reference frame with origin in the molecular center of mass. As 
for shear viscosity, all three independent heat fiow directions , J*', and J|, 
can be used to improve the statistics of Eq. (JU]). 

All results of this study were obtained and are presented in the reduced form, 
i.e. in relation to the molecular parameters size, energy, and mass. The reduced 
transport properties are defined by: D* = D/a^m/e, rj* = rjcr^/^/me, \* = 



Xa^/kB\/m/e. Relevant static thermodynamic properties, temperature T* 
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r/c^/e and number density p* = po"^, are also reduced in the same sense. For 
the sake of brevity, "reduced" will be omitted in the following. 

4 INVESTIGATED MODELS AND STATES 

In the present work, 30 different model fluids were studied, where each fluid 
is fully determined by one combination of elongation L* and quadrupole mo- 
mentum Q*"^. Simulations at 16 liquid state points were carried out for each 
model fluid. In Figs.[T]and[2]four selected systems are shown; they illustrate the 
covered thermodynamic states and the influence of elongation and quadrupole 
momentum on the thermodynamic behavior of the fluids [18j. As Fig. [T]shows, 
critical temperature and density increase with increasing quadrupole momen- 
tum. On the other hand. Fig. [2] shows that these critical properties decrease 
with increasing elongation. Such a behavior can also be seen for linear Kihara 
fluids [31]. 

The studied model fluids have elongations that vary from L* = 0, i.e. spherical 
molecules, to 0.8, i.e. strongly elongated dumbbell-shaped molecules, in six 
steps. The odd value L* = 0.505 was chosen to cover the same model fluids as 
the previous work on vapor-liquid equilibria Five quadrupole momenta 
were studied that range from Q*"^ = to 4 with increments of unity. The upper 
limit of 4 is sufficient to describe strongly quadrupolar real fluids, eg. CO2 with 
Q = -3.7938 DA (Q*2 = 3.3037), C2H2 with Q = 5.0730 DA {Q*^ = 4), or 
C2F4 with Q = -7.0332 DA {Q*^ = 3.9272) [SU]. 

For the sake of consistency, transport properties for spherical fluids (L*=0) 
were treated as two superimposed Lennard- Jones sites. This implies, that the 
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temperature has to be divided by 4 as well as the quadrupole momentum, 
if a direct comparison with a one-center Lennard- Jones fluid is to be made. 
The corresponding conversion of D* , t]*, and A* for these spherical fluids is 
obtained dividing the present data by 2. 

As temperatures and number densities in vapor-liquid equilibrium vary strongly 
with the molecular parameters, it is useful to introduce another reduced form 
for representing the temperature Tji=T* /T* and the density Pr=P* / p*^- Here, 
T* is the critical temperature and p* the critical density of the individual 
2CLJQ fluid; values for T* and p* were taken from |18j. For each fluid, the 
considered reduced temperatures along the bubble line range from Tr=0.6 to 
0.9 with increments of ATr=0.1. In addition to those four points, another 
12 points in the homogeneous liquid region were simulated, cf. Figs. [T] and [21 
These points were selected on isochores starting from each bubble point with 
temperature increments of ATr=0.1. In this way, also isothermal data was 
generated. 

5 SIMULATION DETAILS 

Equilibrium molecular dynamics simulations were performed in a cubic box 
of volume V containing = 500 molecules. The cut-off radius was set to 
Tc = 5cr, otherwise to half of the box length. The molecules in the fluid were 
assumed to have no preferential relative orientations outside of the cut-off 
sphere. For the calculation of the LJ long range corrections, orientational 
averaging was done with equally weighted relative orientations as proposed 
by Lustig [35] . The assumption of no preferential relative orientations beyond 
the cut-off sphere implies that no long range corrections for the quadrupolar 
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interactions are needed since they disappear. This is a reasonable assumption 
as demonstrated by Streett and Tildesley p6j. The simulations were started 
with the molecules in a face centered cubic lattice with random velocities, 
the total momentum of the system was set to zero, and Newton's equations 
of motion were solved with the Gear predictor-corrector integration scheme 
of fifth order [HTj. The time step for this algorithm was set to At-y/e/m/a = 
0.0005. Self-diffusion coefficient, shear viscosity, and thermal conductivity were 
calculated in the microcanonical NVE ensemble using Eqs. (E]) to (JTj). The 
simulations were equilibrated in the canonical NVT ensemble between 100 000 
to 150 000 time steps. After equilibration, the thermostat was turned off and 
the simulation continued in the NVE ensemble where the transport properties 
were calculated. Because of the lack of a thermostat, the temperature was 
fluctuating with a maximum drift of 3 %. 

Statistical uncertainties were estimated using the standard deviation of four 
independent simulations of 3 000 independent autocorrelation functions. In 
order to achieve independence between autocorrelation functions, a time span 
of 0.1 in reduced units was left between consecutive autocorrelation functions. 
This time span was consistent with a decay to less than 1/e of the normalized 
velocity autocorrelation function in several pilot runs. It is a conservative 
choice, when a compromise between simulation time and accuracy has to be 
done. Fig.[3]shows the normalized autocorrelation functions (ACF) from top to 
bottom for self- diffusion coefficient, shear viscosity, and thermal conductivity 
for the 2CLJQ fluid with L*=0.2 and Q*^=l for the bubble point at 70 % of 
its critical temperature {T* = 4.388), where the bubble density is p* = 0.6573. 
The vertical line denotes t*=0.1, the horizontal line the value 1/e. All three 
ACF fulfill the criterion of independence between consecutive correlations. 
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although they show a quite different decay. Another important issue is the 
significant length of the autocorrelation functions, or equivalently, how long 
they must be integrated. Fig. H] shows the integrals from top to bottom for self- 
diffusion coefficient, shear viscosity, and thermal conductivity. In principle, the 
integration of the ACF must be carried out until the integral shows stationary 
behavior. In practice, the convergence of the self-diffusion integral is difficult to 
guarantee for all conditions, because of the long time behavior [4J with a decay 
proportional to t^^^. A similar problem is present regarding the autocorrelation 
function for shear viscosity close to the fluid-solid transition pi38ll39f4Uj . Here 
it was handled by a long evaluation {At* = 1.25) of the ACF, as can be seen 
in Fig. m The integrals converge to the final value at around 0.5, afterwards 
the ACF fluctuate around zero, cf. Fig. [3l without effective contribution to 
the integrals. 

6 RESULTS 

In this section, the simulation results for the transport coefficients are pre- 
sented. Numerical data for self-diffusion coefficient, shear viscosity, and ther- 
mal conductivity are given in Table [T] for elongations from L*=0 to 0.8 and 
quadrupole momenta from Q*'^=0 to 4. All data in Table [T] correspond to state 
points along the bubble line for reduced temperatures of T/j=0.6, 0.7, 0.8, and 
0.9. The complete data set, with 12 additional state points in the liquid region 
for each fluid, is available in [U] and partially included in Figs. U\ El \lDi [HI \IM 
and [HI The effects of elongation, quadrupole momentum, temperature, and 
density are discussed in the following for each transport coefficient separately. 

The accuracy of the calculated transport properties decreases in the sequence 
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self- diffusion coefficient, shear viscosity, thermal conductivity. The high accu- 
racy of the self- diffusion coefficient, with error bars lower than 3 %, is due 
to its individual nature p]. Shear viscosity and thermal conductivity are col- 
lective properties, consequently they show for the same simulation time and 
system size larger uncertainties, that are around 8 and 12 %, respectively. In 
most simulations of the present work the autocorrelation functions of thermal 
conductivity decay faster than those for shear viscosity, but fluctuate more. 
Figs. [3] and H] illustrate this. 

Other factors that influence the accuracy of the reported data are elongation 
and quadrupole momentum. In particular, at low temperatures, for fluids with 
large anisotropy and strong quadrupole momentum, the transport coefficients 
show larger simulation uncertainties. 

In the following, the results are discussed for nine selected fluids, covering the 
whole range of the two molecular parameters, from spherical (L*=0) over elon- 
gated (L*=0.505) to strongly elongated (L*=0.8) fluids with varying quadru- 
pole momentum of Q*^=0, 2, and 4. A subset of six fluids is taken in some 
cases only due to graphical reasons. 

6.1 Self- diffusion coefficient 

Figs. [5] and [6] illustrate the self-diffusion coefficient along the bubble line for 
nine selected fluids. The results can either be discussed in terms of reduced 
density pn as in Fig. or in terms of number density p* as in Fig. [HI From Fig. 
|5] it can be seen that the regarded range of reduced density is similar for all 
fluids, but signiflcant deviations from the principle of corresponding states are 
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present also for the density. At constant Tr, it can be discerned that the self- 
diffusion coefficient decreases with both increasing elongation and quadrupole 
momentum. A better visibility of the data (which is even more needed for 
the less accurate properties shear viscosity and thermal conductivity) is ob- 
tained when plotted over number density in Fig. [HI Therefore, this graphical 
representation is preferred in the following. 

As Fig. [6] shows, D* decreases with increasing number density along the bub- 
ble line (where with increasing density also the temperature decreases). It 
is an important result of the present study that the self-diffusion coefficient 
lies roughly along the same line for a given elongation, independent of the 
quadrupole momentum. 

Fig. [7] shows the dependence of D* on number density in the homogeneous liq- 
uid region at a constant reduced temperature of Tr=0.9. Note that the density 
range is the same as in Fig.O Along this isotherm D* decreases slightly hyper- 
bolic with increasing density, resembling the behavior of -D* along bubble lines 
for a given elongation. Comparing D* along bubble lines with isothermal data 
for an identical density variation, it is found that the density effect dominates 
with a contribution of 80 %. 

Fig. [H] shows the dependence of the self-diffusion coefficient on temperature at 
different constant densities for a subset of six fluids. The isochores correspond 
to bubble densities at the reduced temperature Tr=0.6, cf. Figs. [1] and O 
which have similar values in terms of p^. Along the isochores, the self-diffusion 
coefficient increases linearly with increasing temperature. The gradients with 
respect to reduced temperature are almost constant for a given elongation, 
where the slope is less steep for more elongated fluids. Such a linear dependence 
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of D* on temperature has also been reported by other authors for Lennard- 
Jones fluids [42j , Kihara fluids [7] , and two-center Lennard- Jones fluids [13] . 



6.2 Shear viscosity 

Fig. [9] illustrates the shear viscosity along the bubble line for the nine selected 
fluids. At constant Tr, it is found that the shear viscosity decreases with 
increasing elongation but increases with increasing quadrupole momentum. 
Again it is found that the results for a given elongation lie roughly along one 
line independent of Q*^, where, as expected, they increase with increasing 
number density. 

The density dependence of shear viscosity in the homogeneous liquid region is 
illustrated in Fig. [10] at Tr=0.9. Comparing the variation of r]* along bubble 
lines and along isotherms in the same way as for D*, it is found for non-polar 
fluids that the density effect is responsible for about 80 % of the increase of 
T]* along the bubble line. For quadrupolar fluids, however, the temperature 
influence becomes more important and its contribution is about 40 %. 

Fig. [m shows the dependence of shear viscosity on reduced temperature for 
a subset of six fluids along isochores with similar values in terms of p/j. 
As expected, shear viscosity decreases with increasing temperature. Strongly 
quadrupolar fluids, with an about threefold higher shear viscosity in the cold 
liquid, are more sensitive to temperature, exhibiting larger gradients. 
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6.3 Thermal conductivity 



Fig. [12] illustrates the thermal conductivity along the bubble line. Again, the 
data lie roughly along single lines for a given elongation, but considering sim- 
ulation uncertainties not more than a linear dependence can be discerned. 
Thermal conductivity has the same basic trends like shear viscosity as it de- 
creases with increasing elongation but increases with increasing quadrupole 
momentum at constant Tr. 

Fig. dni shows the density dependence in the homogeneous liquid at Tr=0.9. 
It can be seen that the curves resemble those along bubble lines, underlining 
the dominating effect of density there. Similar results have been reported by 
Tokumasu et al. |44J who studied the non-polar 2CL J potential but at different 
thermodynamic conditions. In their analysis, Tokumasu et al. reduced A* by 
critical temperature and critical density, to isolate the effect of elongation and 
found that this type of reduced thermal conductivity increases with increasing 
elongation. 

Fig. [TH shows isochoric data with similar values in terms of pr for a subset of 
six fluids, where the effect of temperature on A* is small. Taking the statistical 
uncertainty and the scatter into account, hardly any trend can be discerned. 
Experimental results [^HIHB] show that thermal conductivity at constant den- 
sity increases with increasing temperature, but the variation is very small in 
the liquid region. Moreover, the increase of A* with increasing quadrupole 
momentum can be seen. 
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7 CONCLUSION 



Equilibrium molecular dynamics simulation and the Green- Kubo formalism 
were used to calculate self-diffusion coefficient, shear viscosity, and thermal 
conductivity for 30 different anisotropic and quadrupolar model fluids. A com- 
prehensive data set was obtained for each fluid and property that covers a 
substantial part of the liquid state. The statistical uncertainty of the reported 
data varies according to transport property. For self-diffusion coefficient data, 
it is less than 3 %, for shear viscosity and thermal conductivity it is around 8 
and 12 %, respectively. 

The three transport properties are dominated in the investigated liquid region 
by the density: saturated liquid and isothermal data for ffuids with a given 
elongation but varying quadrupole momentum lie roughly along single lines 
when plotted over number density. 

However, all transport properties on the bubble line at a constant reduced tem- 
perature are lower for ffuids with larger elongation. An increasing quadrupole 
momentum also leads to a lower self-diffusion coefficient, the opposite is found 
for shear viscosity and thermal conductivity. 

Temperature inffuences all transport properties less than density. As expected, 
along isochores, the self-diffusion coefficient increases with temperature, the 
shear viscosity decreases, and for the thermal conductivity hardly any varia- 
tion can be discerned. 
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List of symbols 

a interaction site index 

b interaction site index 

c short notation for a trigonometric function 

D self-diffusion coefficient 

i molecule index 

j molecule index 

Jp element of the microscopic stress tensor 

Jg element of the microscopic heat flow 

k interaction site index 

k molecule index 

ks Boltzmann constant 

I interaction site index 

L molecular elongation 

m molecular mass 

N number of molecules 

Q molecular quadrupole momentum 

r site-site distance 

Tc center-center cut-off radius 

s short notation for a trigonometric function 

t time 

T temperature 

u pair potential 

V element of the velocity vector 

V molar volume 
A increment 

At integration time step 

e Lennard- Jones energy parameter 

eo Electric constant 
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7] shear viscosity 

9 angle of nutation 

A thermal conductivity 

p density 

a Lennard- Jones size parameter 

(f) azimuthal angle 

Vector properties 



I matrix of angular momentum of inertia 

Jp microscopic stress tensor 

Jg microscopic heat flow vector 

r distance vector 

V velocity vector 

w angular velocity vector 

r torque vector 



Subscript 



a interaction site index 

b interaction site index 

c property at critical point 

i molecule index 
molecule index 

Q point quadrupole 

R property reduced by critical value 



2CLJQ two-center Lennard- Jones plus point quadrupole 

Superscript 

X cartesian direction 
y cartesian direction 
z cartesian direction 

* property reduced by molecular parameters 
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Table I 

Transport coefficients along tlic bubble line for 30 2CLJQ fluids of different elon- 
gation L* and quadrupole momentum Q*'^. The numbers in parentheses denote the 
uncertainty in the last digits. 



L*=0 T* p* 


D* ri* X* 


Q*2^0 3.156 0.8062 
3.681 0.7453 
4.255 0.6735 
4.674 0.5838 


0.095(2) 4.83(9) 13.4(6) 
0.155(4) 3.00(35) 11.0(13) 
0.250(6) 2.00(24) 8.65(77) 
0.388(5) 1.52(16) 6.37(51) 


g*2=i 3.132 0.8201 
3.712 0.7575 
4.187 0.6851 
4.773 0.5934 


0.080(1) 5.31(35) 11.7(12) 
0.143(2) 3.62(23) 11.2(12) 

0.228(4) 2.26(18) 9.8(12) 
0.368(2) 1.58(6) 6.44(78) 


Q*2=2 3.398 0.8483 
3.925 0.7819 
4.538 0.7041 
4.946 0.6055 


0.062(3) 7.67(46) 13.0(19) 
0.119(1) 4.45(34) 10.7(22) 
0.205(2) 2.64(25) 10.32(82) 
0.339(3) 1.79(16) 7.07(11) 


g*2=3 3.505 0.8796 
4.106 0.8099 
4.741 0.7292 

5.341 0.6272 


0.047(1) 9.70(64) 15.4(25) 
0.099(2) 5.83(58) 12.4(13) 
0.179(1) 3.28(62) 10.4(19) 
0.310(4) 1.93(26) 7.9(12) 


Q*2^4 3g4^ 0.9143 
4.476 0.8430 
5.164 0.7609 
5.742 0.6579 


0.038(1) 14.02(86) 18.3(19) 
0.081(1) 6.87(64) 15.9(17) 
0.153(3) 3.85(21) 11.9(14) 
0.275(4) 2.54(27) 9.8(15) 
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Table I 
Continued. 



L*=0.2 T* p* 


D* r)* X* 


g*2=0 2.589 0.7114 
3.015 0.6573 
3.441 0.5946 

3.874 0.5144 


0.089(1) 3.64(16) 11.58(59) 
0.147(1) 2.63(7) 9.69(74) 
0.228(6) 1.74(11) 7.67(83) 

0.353(7) 1.27(5) 5.90(30) 


Q*2=i 2.625 0.7203 
3.070 0.6644 

3.496 0.5998 
3.941 0.5185 


0.083(1) 4.19(21) 12.1(10) 

0.142(2) 2.86(23) 11.7(13) 
0.226(4) 1.89(9) 8.54(57) 
0.350(9) 1.19(7) 6.74(28) 


Q*2=2 2.722 0.7420 
3.195 0.6833 
3.659 0.6167 
4.072 0.5322 


0.075(1) 4.78(36) 17.1(14) 
0.133(2) 2.97(12) 15.4(15) 
0.215(4) 2.12(14) 11.76(93) 
0.338(8) 1.34(8) 7.74(47) 


Q*2=3 2.877 0.7683 
3.393 0.7085 
3.856 0.6397 
4.318 0.5535 


0.067(1) 5.48(21) 21.0(29) 
0.123(1) 3.67(3) 19.8(16) 
0.203(5) 2.40(7) 15.6(18) 
0.326(6) 1.53(16) 9.51(10) 


Q*2=4 3.054 0.7939 
3.642 0.7311 
4.103 0.6596 
4.600 0.5686 


0.026(3) 11.21(42) 23.3(31) 
0.117(2) 4.06(14) 22.5(24) 
0.196(2) 2.44(8) 18.4(21) 
0.324(10) 1.71(10) 12.0(21) 
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Table I 
Continued. 



L*=OA T* p* 


D* T]* A* 


Q*2=o 1.893 0.5808 
2.232 0.5365 
2.536 0.4853 

2.858 0.-1185 


0.094(2) 2.51(20) 9.36(45) 
0.147(1) 1.82(9) 8.90(36) 
0.219(5) 1.32(9) 6.67(57) 

0.32G(3) 0.89(5) 4.92(47) 


Q*^=l 1.925 0.5879 
2.239 0.5426 
2.573 0.4913 
2.869 0.4252 


0.084(1) 2.82(17) 10.12(60) 
0.134(1) 2.03(7) 8.97(38) 
0.206(1) 1.35(5) 7.39(55) 
0.311(2) 0.92(8) 5.24(84) 


g*2=2 1.999 0.6025 
2.318 0.5555 
2.663 0.5008 
2.987 0.431 


0.072(1) 3.36(11) 11.11(78) 
0.121(1) 2.23(6) 9.58(33) 
0.194(2) 1.53(8) 7.60(55) 
0.304(4) 0.97(5) 6.02(37) 


Q*2=3 2.096 0.6209 
2.447 0.5717 
2.794 0.5154 
3.134 0.4428 


0.061(1) 3.98(14) 14.9(11) 
0.109(1) 2.64(3) 10.52(82) 
0.180(1) 1.71(7) 8.52(38) 
0.290(5) 1.07(6) 6.37(56) 


Q*2=4 2.213 0.6396 
2.599 0.5884 
2.972 0.5307 
3.350 0.4554 


0.051(6) 2.72(13) 14.04(71) 
0.098(2) 2.95(7) 13.3(12) 
0.168(1) 1.90(13) 10.51(17) 
0.281(3) 1.21(4) 7.61(41) 
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Table I 
Continued. 



L*=0.505 T* p* 


D* ri* X* 


g*2=o 1.638 0.5291 
1.913 0.4891 
2.190 0.4431 

2.476 0.3835 


0.095(1) 2.02(6) 9.13(84) 
0.141(2) 1.55(8) 7.62(42) 
0.206(1) 1.10(6) 5.9(10) 

0.305(1) 0.79(4) 4.43(18) 


Q*^=l 1.652 0.5349 
1.924 0.4942 
2.187 0.4474 
2.509 0.3873 


0.083(1) 2.29(14) 9.97(74) 
0.130(1) 1.67(6) 8.08(68) 
0.193(1) 1.19(4) 6.48(18) 
0.295(3) 0.81(5) 4.83(36) 


Q*2=2 1.728 0.5476 
2.029 0.5049 
2.288 0.4548 
2.584 0.3926 


0.071(2) 2.80(10) 10.8(11) 
0.118(1) 1.96(18) 8.70(92) 
0.183(3) 1.34(13) 7.04(71) 
0.283(1) 0.85(4) 4.93(40) 


Q*2=3 1.813 0.5643 
2.102 0.5193 
2.393 0.4691 
2.692 0.4041 


0.058(1) 3.67(17) 11.5(14) 
0.102(1) 2.40(13) 10.6(11) 
0.166(2) 1.52(11) 7.49(65) 
0.266(4) 0.95(10) 5.52(34) 


Q*2^4 ^ 922 0.5803 
2.252 0.5335 
2.541 0.4809 
2.885 0.4135 


0.043(1) 5.28(74) 14.94(28) 
0.092(2) 2.66(13) 11.39(60) 
0.155(1) 1.71(10) 8.52(58) 
0.259(3) 1.02(4) 5.80(24) 
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Table I 
Continued. 



L*=0.6 T* p* 


D* T]* X* 


Q*2=o 1.475 0.4900 
1.726 0.4521 
1.948 0.4088 

2.211 0.3520 


0.094(1) 1.76(8) 9.39(96) 
0.140(1) 1.37(8) 7.67(37) 
0.199(2) 1.02(6) 6.00(68) 

0.294(3) 0.6G(G) 4.13(45) 


g*2=i 1.490 0.4947 
1.731 0.4563 
2.011 0.4116 
2.233 0.3519 


0.084(1) 2.01(8) 8.92(39) 
0.128(1) 1.45(8) 7.46(59) 
0.193(2) 1.10(4) 6.09(37) 
0.287(1) 0.70(3) 3.85(27) 


g*2=2 1.552 0.5083 
1.810 0.4682 
2.080 0.4214 
2.314 0.3622 


0.069(1) 2.42(11) 10.06(31) 
0.112(1) 1.72(9) 8.07(74) 
0.177(2) 1.19(14) 6.09(35) 
0.268(1) 0.79(4) 3.6(11) 


Q*2=3 1.610 0.5239 
1.879 0.4819 
2.137 0.4349 
2.435 0.3758 


0.055(1) 3.04(26) 11.25(61) 
0.096(1) 2.05(6) 9.40(67) 
0.156(1) 1.38(9) 6.95(67) 
0.250(3) 0.90(4) 4.97(38) 


Q*2=4 1.725 0.5381 
2.023 0.4944 
2.284 0.4452 
2.608 0.3812 


0.040(1) 4.84(23) 10.0(24) 
0.087(1) 2.37(6) 10.23(71) 
0.147(2) 1.52(5) 7.85(70) 
0.248(1) 0.99(7) 5.72(25) 
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Table I 
Continued. 



L*=0.8 T* p* 


D* T]* A* 


Q*2=o 1.234 0.4302 
1.426 0.3956 
1.632 0.3568 

1.856 0.3051 


0.085(1) 1.51(5) 7.38(51) 
0.126(1) 1.13(4) 6.21(28) 
0.181(1) 0.85(2) 4.93(13) 

0.267(3) 0.58(5) 3.64(23) 


g*2=i 1.246 0.4364 
1.458 0.4016 
1.686 0.3602 
1.885 0.3059 


0.074(1) 1.69(1) 8.09(50) 
0.115(1) 1.23(7) 6.86(51) 
0.173(1) 0.89(8) 5.56(41) 
0.259(1) 0.60(4) 3.64(30) 


g*2=2 1.308 0.4513 
1.514 0.4143 
1.712 0.3734 
1.933 0.3207 


0.060(1) 2.23(6) 9.52(13) 
0.099(2) 1.46(8) 7.22(22) 
0.150(2) 1.02(3) 5.74(9) 
0.233(1) 0.70(2) 4.04(24) 


Q*2=3 1.352 0.4666 
1.579 0.4290 
1.801 0.3858 
2.045 0.3319 


0.047(1) 2.91(12) 9.02(12) 
0.082(3) 1.85(5) 7.42(29) 
0.136(2) 1.24(6) 6.19(40) 
0.219(1) 0.78(6) 4.42(35) 


Q*2=4 1.447 0.4800 
1.674 0.4416 
1.914 0.3975 
2.187 0.3402 


0.032(1) 4.92(47) 10.81(15) 
0.073(1) 2.18(7) 8.87(40) 
0.125(2) 1.42(1) 7.22(39) 
0.212(1) 0.90(5) 4.65(41) 
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List of Figures 



1 Phase diagrams for two selected elongated 2CLJQ fluids 
(L*=0.2) where one is non-polar and the other strongly 
quadrupolar. Saturated densities, taken from [19], are 
represented by the lines joining at the critical point depicted 
by •. The investigated state points are indicated by o for 
Q*2=o and by A for Q*^=A. 

2 Phase diagrams for two selected quadrupolar 2CLJQ fluids 
((5*^=1) where one is spherical and the other strongly 
elongated. Saturated densities, taken from [19], are represented 
by the lines joining at the critical point indicated by •. The 
investigated state points are indicated by o for L*=0 and by 
A for L*=0.8. 

3 Autocorrelation functions (ACF) for self-diffusion coefficient, 
shear viscosity, and thermal conductivity for the 2CLJQ 
ffuid with L*=0.2 and Q*^=l at T*=3.0295 and p*=0.6573. 
The vertical line denotes t*=0.1. At this time all normalized 
autocorrelation functions have decayed to much less than 1/e 
denoted by the horizontal lines. 

4 Normalized integrals of autocorrelation functions, i.e. 
self-diffusion coefficient, shear viscosity, and thermal 
conductivity. 
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5 Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=0.8, full symbols) 2CLJQ fluids over reduced density 
along bubble lines. Reduced temperatures vary from rR=0.6 
to 0.9. Lines are guides for the eye. 

6 Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=0.8, full symbols) 2CLJQ fluids over number density 
along bubble lines. Reduced temperatures vary from Tr=0.6 
to 0.9. Lines are guides for the eye. 

7 Self-diffusion coefficient of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=0.8, full symbols) 2CLJQ fluids over number density in 
the homogeneous liquid at Tr=0.9. Lines are guides for the 
eye. 

8 Self-diffusion coefficient of spherical (L*=0, empty symbols) 
and strongly elongated (L*=0.8, full symbols) 2CLJQ fluids 
over reduced temperature in the homogeneous liquid along 
different isochores. o: p*=0.8062. A: p*=0.8483, 0: p*=0.9143, 
• : p*=0.4302. A: p*=0.4513, ♦: p*=0.4800. Lines are guides for 
the eye. 
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Shear viscosity of spherical (L*=0, empty symbols), elongated 

(L*=0.505, grey symbols), and strongly elongated (L*=0.8, 

full symbols) 2CLJQ fluids over number density along bubble 

lines. Reduced temperatures vary from Tr=0.6 to 0.9. Lines 

are guides for the eye. |40 



10 Shear viscosity of spherical (L*=0, empty symbols), elongated 
(L*=0.505, grey symbols), and strongly elongated (L*=0.8, 
full symbols) 2CLJQ fluids over number density in the 
homogeneous liquid at Tr=0.9. Lines are guides for the eye. |41 



11 Shear viscosity of spherical (L*=0, empty symbols) and 

strongly elongated (L*=0.8, full symbols) 2CLJQ fluids over 
reduced temperature in the homogeneous liquid along different 
isochores. o: p*=0.8062. A: p*=0.8483, 0: p*=0.9143, •: 
p*=0.4302. A: p*=0.4513, ♦: p*=0.4800. Lines are guides for 
the eye. |42 



12 Thermal conductivity of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=0.8, full symbols) 2CLJQ fluids over number density 
along bubble lines. Reduced temperatures vary from Tr=0.6 
to 0.9. Lines are guides for the eye. |43 
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13 Thermal conductivity of spherical (L*=0, empty symbols), 
elongated (L*=0.505, grey symbols), and strongly elongated 
(L*=0.8, full symbols) 2CLJQ fluids over number density in 
the homogeneous liquid at Tr=0.9. Lines are guides for the 

eye. M 

14 Thermal conductivity of spherical (L*=0, empty symbols) 
and strongly elongated (L*=0.8, full symbols) 2CLJQ fluids 
over reduced temperature in the homogeneous liquid along 
different isochores. o: p*=0.8062. A: p*=0.8483, 0: p*=0.9143, 
• : p*=0.4302. A: p*=0.4513, ♦: p*=0.4800. Lines are guides for 

the eye. |45 
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Fig. 1. 
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Fig. 2. 
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Fig. 3. 
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Fig. 4. 
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Fig. 5. 
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Fig. 6. 




37 



Fig. 7. 
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Fig. 8. 
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Fig. 9. 
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Fig. 10. 
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Fig. 11. 
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Fig. 12. 
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Fig. 13. 
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Fig. 14. 
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